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We study in-gap electronic states induced by a nonmagnetic defect with short-range potential in 
two-dimensional topological insulators and trace their evolution as the distance between the defect 
and the boundary changes. The defect located far from the boundary is found to produce two 
bound states independently of the sign of its potential. The states are classified as electronlike and 
holelike. Each of these states can have two types of the spatial distribution of the electron density. 

The first-type states have a maximum of the density in the center and the second-type ones have 
a minimum. When the defect is coupled with the boundary the bound states are transformed 
correspondingly into resonances of two types and take up the form of the edge states flowing around 
the defect. Under certain conditions, two resonances interfere giving rise to the formation of a 
bound state embedded into the continuum spectrum of the edge states flowing around the defect. 

We calculate the spatial distribution of the electron density in the edge states flowing around the 
defect and estimate the charge accumulated near the defect. The current density field of the edge 
states flowing around the defect contains two components one of which flows around the defect and 
the other circulates around it. 


I. INTRODUCTION 

The presence of gapless edge states at the interface of 
topologically non-equivalent crystals is a hallmark of two- 
dimensional (2D) topological insulators (TIsj- 1 . In these 
states the electrons move along the boundary and their 
spin is locked to the momentum because of strong spin- 
orbit interaction. Such helical edge states are protected 
against scattering by weak non-magnetic impurities and 
disorders. Nevertheless, experiments reveal a noticeable 
backscatterin g of e lectrons 2 the mechanism of which is 

not yet knowipMI. 

Backscattering of electrons in the edge states can oc¬ 
cur as a result of an inelastic process due to electron- 
electron interactions and the presence of a defect po¬ 
tential 11 . The effect of the electron-electron interaction 
in the vicinity of the defect essentially depends on the 
charge and spin structure of the electron cloud which 
forms near it. In this regard, of great importance is the 
question about the electronic states induced by impuri¬ 
ties and other structural imperfections, especially in the 
case where the defect is located near the boundary. One¬ 
dimensional models of coupling between the edge states 
and the defect turn out to be insufficient to describe the 
experiments^®!!. 

Electronic states induced by a defect were studied for 
three dimensional (3D) TIs where the defect is located 
on the surface. In this case, the electron cloud around 
the defect is formed by 2D electronic states propagating 
along the surface. Their interference leads to a variety of 
the electron density configurations— 17 and even to the 
changes in the surface state spectrrP®®. 

In 2D TIs, the electron cloud around a defect also exists 
but its structure is substantially different since the elec¬ 
tron density configuration is formed mainly by evanescent 
modes decaying in the plane. It is essential that the elec¬ 
tron cloud can not be described within a one-dimensional 


(ID) model. Electronic structures formed in this case are 
currently poorly understood. 

Defect-induced electronic states in 2D TIs were stud¬ 
ied mostly in the case where the defect is located deep in 
the bulk and decoupled from the boundary. In Ref. [23 
the defect was considered as a hole, at the edges of which 
the wave function is zero. In this case, the bound states 
are in essence the edge states circulating around the hole 
with quantized angular momentum. Although this model 
captures some properties of the defect-induced states, it 
is far from reality. Under realistic conditions the wave 
function is not zero in the defect. The bound states ap¬ 
pearing in the Gaussian potential were investi gated nu¬ 
merically for a number of material parameter j 21 l 22 l , but 
no general conclusions were made about their spectra, 
the electronic structure, and the conditions under which 
they exist. 

Defects interacting with the boundaries were studied 
in the case of a slab of 2D TI. In this case, the defect 
is coupled with two boundaries. Numerical calculations 
with using Green’s function method combined with tight- 
binding approactPHl have shown that the bound state 
spectrum differs from that in the continuous model. Par¬ 
ticularly, it contains two states bound on one defect with 
short range potential rather than one state as in the con¬ 
tinuous model 21 . 

In recent work^i we investigated analytically the 
bound states induced by a non-magnetic defect in the 
bulk of 2D TIs for defects with short-range potential. 
It turned out that the defect creates two bound states 
which are classified as electronlike and holelike. This is 
in contrast to the defects in topologically trivial insula¬ 
tors where only one bound state exists in a short-range 
potential. The bound states exist for both positive and 
negative potentials. In turn these states can be also of 
two types depending on whether the electron density has 
a maximum or a minimum in the point of the defect lo- 
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cation. Another interesting feature of the 2D TIs is an 
unusual dependence of the bound-state energies on the 
defect potential. As the potential increases, the energies 
of both electronlike and holelike states tend correspond¬ 
ingly to two different limiting values, which lie within the 
gap. 

In this paper we address the general problem of a de¬ 
fect coupled with the edge states in 2D TIs. We clarify 
how the bulk bound states are modified with decreasing 
the distance between the defect and the boundary and 
how the edge states are distorted by the defect. It turns 
out that the edge states and the bulk bound states trans¬ 
form into a set of eigenstates which have the form of the 
edge states flowing around the defect. These states have 
resonances of the electron density in the vicinity of the 
defect when the energy is close to the energy of the bulk 
bound states. Correspondingly, there are two types of 
the resonances. 

Under certain conditions two resonances of different 
types can interfere with each other giving rise to the for¬ 
mation of a bound state with localized wave function in 
the continuum of the edge states. 

We study the spatial distribution of the electron den¬ 
sity and current density in the states flowing around the 
defect and estimate the charge accumulated near the de¬ 
fect at a given Fermi energy. 

The outline of the paper is as follows. In Section [II] 
we present analytical calculations showing the presence 
of the states flowing around the defect and the bound 
states in the continuum. Section |III| gives the detailed 
results for the bound states in the bulk of 2D TI. Sec¬ 
tion IIVI deals with the electron states in the case where 
the defect is located at a finite distance from the bound¬ 
ary. In Sec.[Vj we study the electron density distributions 
for resonant states, estimate the excess electron density 
accumulated near the defect and consider the patterns of 
the electron current. We finish the paper with a discus¬ 
sion and conclusions in Sec. eh 


II. BOUND STATES AND EDGE STATES 
FLOWING AROUND THE DEFECT 


Our study is based on the model of the 2D TIs 
proposed by Bernevig, Hughes, and Zhang (BHZ) for 
HgTe/CdTe quantum wells^. The 2D TI is described 
by the Hamiltonian 


H 0 = 



° 1 


( 1 ) 


where k is momentum operator and 


fM — {B+D)k 2 A(k x +ik y ) \ 

y A(k x — iky) —M+(B — D)k 2 J ’ 


( 2 ) 



FIG. 1. (Color online) A schematic view of a defect located at 
the distance yo from the boundary of the 2D TI. The darkened 
area shows the particle cloud near the defect. Lines represent 
the particle flows. 


M, B,D < 0, and A > 0. The basis set of wave functions 
is {\Ei t), l-Hi t ),\Ei I ),\Hi 4-)} where \E X f) and \E 1 J,) 
are superpositions of the electron states of s-type and 
light-hole states of p-type with spin up and spin down; 
l-ffi f) and |Hi 4-) are the heavy-hole p-type states with 
opposite spins. In what follows we will restrict ourselves 
by considering the symmetric model where D = 0. 

Let us use the Cartesian coordinates, with the x axis 
coinciding with the boundary (Fig. [I]). The TI lies at 
y > 0 and the defect is located in the point x = 0,p = 
yo . We consider the defect described by a potential 
V(x,y — yo) localized in a small region. Since the defect 
is non-magnetic, the total Hamiltonian H 0 + V(x, y~yo) 
is separated into spin blocks. For spin-up electrons, the 
Schrodinger equation reads as 

[Eao - h(k)] ty{x,y) = a 0 V(x,y- y 0 )^(x,y ), (3) 

where do is a 2 x 2 unit matrix, \&(a;,y) is a spinor 
(il>i(x,y),ijj 2 (x,y)) . The wave functions are supposed 
to vanish at y —X oo and equal zero at y = 0. 

In what follows we will use dimensionless variables 

e = E/\M\, {x',y'} = {x,y}^/M/ B , a = A/y/MB , , . 

v(x',y') = V(x,y)/\B\, b = yo^/M/B, 

and for convenience will omit the prime in the variables 

x',y'. 

The 2D problem ©> can be solved by using the Fourier 
and Laplace transforms over x and y: 

/ OO nOO 

dxe~ lkx / dye~ PVi S>(x, y) , (5) 

-oo J 0 

OO C+200 

*(*,V) = J ^ kx J ( 6 ) 

— OO c—200 

When applying this transformation to Eq. |3| one needs 
to calculate the Fourier and Laplace transforms of the 
product v(x, y)^(x, y). We suppose that the region, 
where the defect potential is localized, is small compared 
with the characteristic length scale of the wave function. 
In this case the integral can be approximated as 


with M, A, B and D being the model param eters. 
The topological phase is realized when MB > () ! 25 l 26 l 
In the case of the HgTe/CdTe wells, the parameters 


OO OO 

J dxe~ lkx J dye~ pv v(x,y)'i/(x,y) 

— oo 0 


v(k,p)e bp tH, (7) 
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where \P = = 0,y = b) is the wave function at 

the defect position and v(k,p) is the Fourier and Laplace 
transforms of v(x,y). In such a way we arrive at the 
following equation: 

[e - h(k,p)]^(k,p) = a z $(k) + a 0 v{k,p)e~ bp 'S> . (8) 

Here, [e — h(k,p)\ is the matrix with elements (e, fc, p): 

an = £ + 1 - k 2 + p 2 , a 12 = -a(k +p), 
a 2 i = -a(k-p), a 2 2 = £ - 1 + k 2 -p 2 , 

<r z is the Pauli matrix, 4>(fc) is the Fourier transform of 
the normal derivative of ^(x,y) at the boundary: 

$(fc) = [ dxe~ ikx ° 

J —oo uy 



equation A(e, fc,pi i2 ) = 0. In such a way we arrive at the 
following equation: 


A(e, k)$(k) + B(e, fc)'P = 0, 


where A(e,k) and B(e,k) are matrices 


A(e, k) 


( a 22 (£, k,pi) a 12 (e,k,pi)\ 
\a 22 (£, k,p 2 ) a 12 (e,k,p 2 )J ’ 


(16) 


(17) 


7?(e, k ) = 

f v(k,pi)a 22 (e, k,p 1 )e~ bpi -vik^^a^e, k,pi)e~ bpi \ 
\y(k,p 2 )a 22 (£,k,p 2 )e~ bp2 ~v(k,p 2 )a 12 (e, k,p 2 )e~ bp2 J ' 

(18) 


Solving Eq. (16) with respect to $(&) we obtain an 
explicit expression for <b(/c): 


We are going to get a system of linear equations for 
the components of the spinor VP which will allow one to 
determine the eigenenergy spectrum. This idea is imple¬ 
mented as follow^ 

Solving Eq. (|8j) with respect to \P(fc, /?) and using 
Eq. we obtain the following expression for the wave 
function: 




CO 



c+ioo 

dk ik f dp e py 

27t 6 J 2nA(e,l,p) 


— OO c—200 

x [D 0 (e, k,p)${k) +v(k,p)e~ bp D 1 (£, k,p)^] , (11) 

where A(e, k,p) is the determinant of the matrix in the 
left-hand side of Eq. Q which has the form: 

A(e, k,p) = — [p 2 - Pi(e, k)] [p 2 - p\(e , k)] , (12) 

with 


Pi j2 (e, k) = \Jk 2 + a 2 /2 — 1 ± \Ja 2 (a 2 — 4)/4 

(13) 

and Repi, 2 (£, k) > 0; D 0 (e,k,p ) and Di(e,k,p) are the 
following matrices: 


a 22 (e,k,p) 

a 12 (e,k,p) \ 

(14) 

—a 2 i(e, k,p) 

-au(e,k,p)J ’ 

a 22 (£, k, p) 

-a 12 (e,k,p)\ 

(15) 

-a 2i (e,k,p) 

an (e,k,p) J ' 


m = _ A '( £ ’Q B (e> k ) v + C(e ) x ( £, k)5[k k 0 (e)], 

Ai(e, k) 

(19) 

where Ai(e, k) is the determinant of the matrix A(s,k) 
and A!(e, k ) is the following matrix: 


( ai 2 (£, k,p 2 ) -a 12 (e, k,p ± )\ 

a 22 (£, k,p 2 ) a 22 (e,k,pi) J 


( 20 ) 


The second term in Eq. ( |19| arises because of the singu¬ 
larity of the matrix A{e, k ) in accordance with the general 
theory of singular matriceiPS. It describes the contribu¬ 
tion of the edge states in the pure TI into the electronic 
states formed in the presence of the defect. k 0 (e) is a 
root of the determinant Ai(e, k) which gives exactly the 
spectrum of the edge states in the absence of the defect: 


ko{e) = -- • 


( 21 ) 


Further in Eq. (19), the coefficient C(s) is a normal¬ 


ization constant, x(e, k) is a spinor which is expressed 
via the matrix elements a^E, fc,p) at p = pi i2 . Using 
Eqs (13) and (21) it is easy to show that %(e, k) coin¬ 
cides with the spinor of the edge states: 


X = 



( 22 ) 


Let us now apply Eq. (11) to calculate T. To this end, 
we set x = 0 and y = b and exclude $(fc) using Eq. (19). 
Finally, we obtain the equation which determines ’P: 


(ao-/C(£))vP = C(£)^(£) X , (23) 


Let us now turn to the requirement that 'P(a;, y) should 
not diverge in the limit y —> oo. Equation shows 
that T (x,y — > oo) -/» oo when the expression in the 
square brackets equals zero at p = Pi j2 (e, k). This gives 
us two equations that relate $(/c) and d'. Correspond¬ 
ingly, there are four equations for their spinor compo¬ 
nents. However, one can show that only two of these 
equations are independent because the matrix elements 
ciijis, k,p) at p = pi and p = p 2 are connected by joint 


where K(e) and JF(e) are the following matrices 


OO 

m= I 


dk 1 

27t 4a e Ai(£, k) 

ico 

f dp v(k,p) 


k)A'(e, k)B(s, k ) 


27tj A(e, k,p ) 


Di (e,k,p) 


, (24) 
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J'(e) = — V 0 (e,k 0 (£)). 


(25) 


Here a e = ^/a 2 (a 2 /4 — 1) + e 2 and V 0 (e,k) denotes the 
matrix 


o-f>Pl 


o~bp2 


V 0 (£, k) = - D 0 (e, k-pi) - D 0 (e, k,-p 2 ). (26) 


Pi 


P 2 


Equation (23) has solutions of two kinds depending on 


the determinant of the matrix (oo — A(e)) 

A*(e) = (1 - A n (e)) (l - A 22 (A)) - A, 2 (e)/C 21 (e) (27) 


First, if A$(e) ^ 0, the root of Eq. (231 is 
C{e) 


*(e) = 


A*(e) 


[cr 0 - A'(e)] J r (e)x, 


(28) 


where 


w/ \ /l — A 2 2(e) Ai2(e) \ /rjQ-, 

- K {C) = ( C 21 (e) 1-AwJ- (29) 


An alternative is the case where 
A*(e) =0. 


(30) 


Let Eq is a root of this equation. When e = e 0 , Eq. (23) 
has a solution if C = 0. This solution reads as 


* (e) = C bs 


1 


(1 — An) / A12 


(31) 


with the constant C bs being determined by the normal¬ 
ization. 

It is worth noting that in the first case, C(e) should 
turn to zero when e tends to c 0 . Otherwise, \h(a:,y) will 
not be normalized. Thus vE'(e) is determined by Eq. (28) 


for any e ^ eq. But if C(e) is exactly zero, the solution 


T is given by Eq. (|31). 

In order to clarify the nature of these solutions we con¬ 
sider the asymptotics of T(a:, y) at x —> ±oo. The asymp¬ 
totic behavior of 4/(a;, y) is easily found from Eqs (11) and 
(19). It has the following form: 


\1/(a? —► oo, y) ~ 


le 


ikx 


8a, 


gAi 

£ dk 


o~P2V 


V 2 


-D 0 (£,k, -p 2 ) 


o-piy 


Pi 


~Dq(e, k, —pi) 


A'(e, k)B(E, fc)4/(e) 


fe=fco(e) 


(32) 


In the case where 4/(e) is determined by Eq. (28), one 
can show that ^{x —> oo ,y) never equals zero and is 
proportional to exp[zfcoa;]. Hence, these states propagate 
along the edge and flow around the defect. We will call 
them the edge states flowing around the defect. They 
have the continuous spectrum defined by Eq. (21), which 


coincides with the spectrum of the edge states without 


defects. The constant C(e) can be found by appropriate 
normalization. 

At a discrete energy e = e 0 the wave function should 
be square integrable and hence the amplitude in its 
asymptotics, given by Eq. (32), should be zero. Using 


the specific expressions for matrices A'(e, k ) and B(e, k), 
given by Eqs. (20) and (18), one can easily show that 


A'(e, k)B{E 1 k)^/ = 0 if 




(33) 


Thus, when 'L(eo) satisfies Eq. (33), a bound state 
can arise in the continuum of edge states. Comparing 
Eq. (31), which defines 4/(e), and Eq. (33) we arrive at 


the following equation for the elements of the A matrix: 


1 — An(e) — Ai 2 (e) — 0. 


(34) 


Importantly, this equation must be satisfied together 
with Eq. (30) that gives the necessary condition for the 


bound state to exist. At this point, one should take into 
account that the elements of the A matrix depend not 
only on the energy e, but also on the defect potential 
v(x,y). Therefore, the system of Eqs. (30) and ( |34j ) de¬ 
termines the energy £ bs of the bound state in the con¬ 
tinuum and the defect potential v bs at which this state 
arises. 

Following, we present the results of specific calculations 
of the bound states and the states flowing around the 
defect. 


III. BOUND STATES IN THE BULK 

We start by considering the limit of b —> oo, which de¬ 
scribes the bound states for a defect located in the bulk. 
When b 


oo, the right-hand side of Eq. (23) goes to 
zero and the nondiagonal components of the A(e) ma¬ 
trix defined in Eq. (24) also vanish. As a result, Eq. (23) 


decouples into two independent homogeneous equations 
for the components of the spinor 4/ = (i/> 1; ■i/i 2 ) T . Cor¬ 
respondingly, there are two kinds of bound states with 
different pseudospin components of the wave function at 
the defect. 

There is a solution in which ij; 1 ^ 0 and z/> 2 = 0. Since 
"01 corresponds to the \E1) component of the basis set of 
wave functions, the states of this kind can be convention¬ 
ally called the electronlike states. In another solution, in 
contrast = 0 and 0 2 ^ 0. We call them the holelike 
states. 

The eigenenergies of the states of both species are de¬ 
termined by Eq. (30). In the limit b —> oo, Eq. (30) de¬ 


couples into two equations. Correspondingly, there are 
two solutions for electron-like and hole-like states: e e 
and Eh- The bound-state energies depend on the defect 
potential v(x,y ) = vf(x,y). Although the particle-hole 
symmetry is broken due to the defect potential, the fol¬ 
lowing symmetry relation holds for the energies of the 
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electron-like and hole-like bound states: 


e e (v) = -e h (-v). 


(35) 


Specific calculations of the bound-state energies and 
the electron density were carried out for the defect po¬ 
tential of two forms: the Gaussian function v(x, y ) = 
uA 2 /7rexp[— A 2 (x 2 + y 2 )] with the characteristic radius 
A -1 , and the v(x,y) = v/tt6 (a" 2 + y 2 ) with regularizing 
cut-off at A when integrating over the wave vector. Both 
cases give similar results. 

Unusual properties of the bound states in 2D TIs be¬ 
come apparent from the dependence of the bound state 
energies on the defect potential amplitude v. They are 
illustrated in Fig. [2ja) . The energies of both electronlike 
and holelike states have two branches with the quite dif¬ 
ferent dependence of the energy on v. To be specific, we 
consider the electronlike states. One branch, £ e \{v), ap¬ 
pears when the potential is attractive for electrons, v < 0. 
As |u| increases, the bound state |el) appears with the 
energy at the top of the gap. Thereafter, its energy goes 
to the bottom of the gap, reaching asymptotically a lim¬ 
iting value e e . We call such states the states of the first 
type. 

When the potential is repulsive for electrons, there is 
another branch e e 2 (v), which represents the bound states 
of the second type, |e2). With increasing v, the energy 
e e 2 changes from the bottom of the gap to the limit¬ 
ing energy £ e . The holelike states \hl) and \h2) behave 
symmetrically with respect to the electron-like states in 
accordance with Eq. (35): £h(i; 2 ){ v ) = ~ £ e(i, 2 )(~ v )- 

A physical difference between the states of the first 
and second types is seen from the spatial distribution of 
the electron density p = and the pseudospin com¬ 
ponents of the density \ipi\ 2 and j^l 2 - Graphs of the 
radial distribution of the total electron density and the 
pseudospin components are shown in Figs. HM-ge). In 
the states of the first type the density p has a maximum 
in the point of the defect location, while in the second- 
type states the density reaches a minimum at the defect. 
Nevertheless, one should note that in the general case 
the maximum of the density in the center is not neces¬ 
sarily the highest maximum in the radial distribution of 
the density. Under certain conditions, another maximum 
may appear at some distance from the center. 

It is remarkable that in the first-type states, the pseu¬ 
dospin component (electronlike or holelike one) which 
reaches a maximum in the center is exactly that for which 
the defect potential is attractive. The opposite situation 
occurs for the states of the second type. The compo¬ 
nent of the spinor which vanishes in the center is that for 
which the potential is repulsive. 

The existence of two states in a short-range potential 
is a feature of the 2D TIs. In the topologically trivial 
case, where MB < 0, the calculations carried out by 
the same method show that there are also electronlike 
and holelike states, but only one state arises in a given 
potential. The electronlike state exists only at v > 0 and 
the holelike state exists at v < 0. Moreover, the states 





FIG. 2. (Color online) A schematic view of (a) the energy of 
the bound states in the bulk of 2D TI as a function of the 
defect potential, (b-e) the radial distribution of the electron 
density p (thick lines) and the densities of the spinor compo¬ 
nents |i/ , i, 2 | 2 for electronlike and holelike bound states of the 
first and second types. 


occur in a finite range of |u|. In both states, the density 
^14' reaches the maximum in the center, i.e. both states 
are the states of the first type in our classification. The 
second-type states are absent. 

These facts allow one to interpret the presence of two 
bound states in a given potential as a result of a simul¬ 
taneous action of two mechanisms of bound state forma¬ 
tion. The first mechanism is universal: the bound states 
can be formed by the potential attracting the quasipar¬ 
ticles of one of the bands. Another mechanism is specific 
for TIs. It is caused by the formation of an edge state 
circulating around the defect similarly to the edge states 
near the boundary. In a certain sense, the defect effec¬ 
tively creates a boundary condition for the wave function. 
This mechanism was discussed in the literatur^ 2022 . 

The existence of two states agrees qualitatively with 
recent numerical calculations with using a tight-binding 
approach combined with the Green’s function method 2 -^, 
but there are essential discrepancies. The contradictions 
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are clearly seen in the spatial distributions of the total 
density p and the pseudospin components of the densi¬ 
ties, as well as in the dependence of the densities on the 
impurity potential. The results we have obtained here 
very well agree with the direct calculations within the 
continuous model of the isolated defect. 

The bound state energies depend also on the parameter 
a, defined through the parameters of the BHZ model by 
Eq. Q. When a < 2 1 / 2 , the energy gap is less than \M\. 
In this case, two bound states with the energy within this 
reduced gap exist for all values of the defect potential and 
the graphs of £( e ,h){ 1 , 2 ) (v) look as in Fig. B a ) . In the 
case where a > 2 1 / 2 , a qualitative difference arises for 
the energy of the first-type states. These states appear 
when |c| exceeds a threshold value. Below the threshold, 
only the second-type states exist. 

The approach we use, does not allow us to investi¬ 
gate the effect of the shape of the defect potential on 
the bound states. One can only trace how the bound 
state energies change with the localization radius when 
it is small. We considered the case where the localization 
length changes together with the potential amplitude so 
that the integral of the potential over the area remains 
constant. It was found that as the localization length de¬ 
creases, the limiting energies e e ,/i shift slowly (logarith¬ 
mically) to the nearest edges of the gap. In the limiting 
case of A —> 00 , the potential shape becomes the 5 func¬ 
tion and the bound states disappear in agreement with 
the theory of singular potentials.^. 

The energy of the bound states does not depend on the 
spin, which means that there are two states with oppo¬ 
site spins. In these states, the electron current circulates 
around the defect in opposite directions, just as in the 
edge states. However, in contrast to them, each state 
can be occupied by one electron with some spin because 
the capture of another electron is hindered by Coulomb 
repulsion. The possibility of capturing another electron 
with opposite spin requires a separate consideration, tak¬ 
ing into account the interaction between electrons. 


IV. DEFECT-INDUCED STATES NEAR THE 
BOUNDARY 


When the defect is coupled with the boundary, the 
electronic states are classified as the edge states flowing 
around the defect and the bound states in the continuum. 
In this section, we consider the resonances of edge states 
flowing around the defect and find out the conditions 
under which the bound state arises in the continuum of 
the edge states. 

The specific calculations are performed for the 5- 
like potential of the defect in the form v(x,y — b) = 
(v/i t) 5 [x 2 + (y — b ) 2 ] with using a cutoff at A in the mo¬ 
mentum space. This simplification allows us to calculate 
analytically the integrals over p in Eqs. (11) and (l24|. 


The subsequent integration over k is done numerically. 

It is natural to expect that when the defect is located 


far from the boundary, the mixing of the bound state 
and the edge states leads to the formation of a resonance 
with the energy close to the bound-state energy. At a 
finite distance b between the defect and the boundary, 
the resonance broadens and the resonant energy deviates 
from the bound-state energy. It is this mixture of the 
states that forms the edge states flowing around the de¬ 
fect. They are exactly described by Eqs pT) ), p9| ), and 

@. 

Our analysis shows that the resonant energy is very 
close to the energy e 0 defined by the roots of Eq. (30). 
In the limit b —>■ 00 , these roots describe the bulk bound 
states. In the case of arbitrary distance b , it becomes 
essential that the roots are functions of three variables: 
the defect potential v, the distance b , and the material 
parameter a. Particularly, the dependence of £ 0 (^, 6 , a) 
on v is much more complicated than in the case of the 
bulk position of the defect. In the following, we con¬ 
sider the dependence of the resonant energy on all three 
quantities. 


When a > 2 1 / 2 , general patterns of the behavior of 
£o{v) with decreasing b are as follows [see Fig. Ha) for 
illustration]. The limiting energies £ ej h shift from their 
positions at b —> 00 to the edges of the gap so that at 
some finite distance b = b m the limiting energy £ e crosses 
the bottom of the gap and £h crosses the top of the gap. 
Since our calculations are not justified for the energy out¬ 
side the gap, we can only conclude that the roots £( e ,h )2 
corresponding to the second-type resonances disappear 
in the gap when b < b m . Nevertheless, the first-type res¬ 
onances continue to exist, but in a finite interval of v, 
which diminishes with decreasing b. 


The minimum distance b m , above which the resonant 
energy of the second-type states lies in the gap, depends 
on the parameter a. Numerical calculations show that 
b m varies slowly between values 2 and 3 when |a| < 1.5 
passing through a minimum at |a| ~ 1. With increasing 
|a| above 1.5, the distance b m grows sharply, reaches 15 
at |a| = 2.5, and continues to grow further. 

A special situation arises in the region near the inter¬ 
section point of the curves £( e ,h)i( v ) and £(h,e) 2( v ) cal¬ 
culated in the limit b — > 00 . At finite b , the coupling of 
the bound states with the boundary results in a mixing of 
the states of the first and second types. This leads finally 
to an avoided crossing of the energy levels £t e ,h)i( v ) anc l 
£(h,e) 2 (v). An example of such an anticrossing is shown 
in Fig. Hb) for 6 = 5. The anticrossing of the levels 
£ (e,h)i( v ) and £ (h,e)2( v ) occurs only when the distance 6 
is large enough. The minimum distance 6 , above which 
the anticrossing occurs, depends upon the parameter a. 
With decreasing a the minimum distance diminishes. 

When a < 2, the wave vectors of the evanescent states, 
Pi, 2(^1 k) defined by Eq. (13) become complex. How¬ 
ever, until the imaginary part of p\ p is small, no new 
effects appear in the behavior of £0 with changing v. A 
qualitatively new behavior appears when a<l and the 
evanescent states strongly oscillate with distance. Their 
interference in the segment between the defect and the 
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FIG. 3. (Color online) (a) Resonant energies £( e ,h)i{v) and 
£(h,e) 2 {v) as functions of the defect potential for a vari¬ 
ety of distances between the defect and the boundary, b = 
2, 2.6, 3.5, 5. The dotted line shows the region magnified in the 
lower part of the figure, (b) Graphs £( e ,h) i(v) and £(h,e) 2 (v) 
near the anticrossing point. Dashed lines show the crossing 
curves £hi(v) and £ e 2 (v) for the bulk states, b = oo. The 
asterisk shows the energy of the bound state £b a = —0.984 
which exists at Vb s = 1.959 when the defect is located at the 
distance b = 5. The calculations were carried out for a = 2 1 ' 2 . 


boundary results in a nonmonotonic dependence of £q on 
v, as it is illustrated in Fig. [4] An oscillatory compo¬ 
nent of the dependence of £ 0 on v is seen to arise when 
the defect approaches the boundary. The effect can be 
so strong that multiple resonant states at a given v can 
exist under certain conditions. 


These effects could be realized in quantum wells 
InAs/GaSb where the parameters of the BHZ model are 
such that a ~ 0.2-0.5 and the imaginary part of pi^ 
greatly exceeds the real part. 

An unexpected effect of the coupling between the de¬ 
fect and the boundary consists in the appearance of the 
bound state in the continuum of the edge states. In 
Sec. |nj we have shown that the energy £& s of this state 
and the defect potential vi, s , at which it appears, are de¬ 
termined by the system of Eqs. (30) and (341. It is not 
difficult to analyze these equations in the case of large 
b where the defect is weakly coupled with the boundary 
and therefore the non-diagonal elements of the K. matrix 
given by Eq. (24) are small. In this case, the non-diagonal 
terms can be treated perturbatively. Finally, we come to 
the following conclusions: (i) the bound state can arise at 
those e and v which are located on the plane (e, v) in the 
vicinity of the intersection point of the curves £r e ,h)i{ v ) 



-12 -10 -8 -6 -4 -2 0 2 4 6 8 10 12 

defect potential, v 

FIG. 4. (Color online) The resonant energies £( e ,h)i(v) and 
£(h,e) 2 (v) as functions of the defect potential for a variety 
of the distance between the defect and the boundary, b = 
2.5,3,10, in the case where a = 0.4. 


and £(h,e) 2 {v ) calculated for b —>• oo. (ii) in this region 
there is only one solution, which exists if b is large enough. 

Thus, the bound state arises at a definite potential Vb s 
when the defect is located at a distance larger than a 
threshold value. The energy of this state £b s and the de¬ 
fect potential vi ls are close to point where two resonances 
are degenerate. The results of specific numerical calcula¬ 
tions for 6 = 5 and a = 2 1 / 2 are shown in Fig. 0b). The 
bound-state energy is indicated by the asterisk. 

This qualitatively new property of the defect-induced 
states is caused by the presence of two types of the reso¬ 
nant states in 2D TIs. The bound state in the continuum 
arises due to the interference of two resonances tuned by 
changing the defect potential so that they can be driven 
into degeneracy. This mechanism is consistent with the 
theory by Friedrich and Wintgen 29 ) 


V. ELECTRONIC STRUCTURE OF STATES 
FLOWING AROUND THE DEFECT 


In this section, we present the results of our calcu¬ 
lations of the spatial distribution of the electron density 
and current density in the resonant states flowing around 
the defect. 

The electron density in the edge state flowing around 
the defect at a given energy, 


Pe{x,y) = ^l{x,y)^ e (x,y), 


(36) 


is calculated using 4/ e defined by Eq. (11 ) where 4>(£, k) 
and 'F(e) are given by Eqs. ( fl9| ) and ([28]). The results 
obtained for the defect located at the distance 6 = 3 from 
the boundary are presented in Fig. [5j The calculations 
were carried out for the (5-like defect potential. The am¬ 
plitude of the potential, v = 3.6, was chosen such that 
one of the resonances £hi was lying deep in the gap, and 
























ter, while in the second-type resonances the density has a 
minimum. In contrast to the bound states in the bulk, no 
component of the spinor T equals exactly zero because 
the resonant states are a mixture of the bulk bound state 
and the edge states. Nevertheless, one of the spinor com¬ 
ponents remains much smaller than the other. 

Now, it is interesting to clarify what charge of electrons 
is accumulated near the defect. The excess electron den¬ 
sity at an energy level e is evaluated as 

Ap e ( x, y) = p e (x, y ) - p s (x, y) ^ =0 , (37) 

where the second term in the right hand side is the den¬ 
sity in the absence of the defect. To simplify the pre¬ 
sentation of the results, we will characterize the excess 
number of electrons by an integral value A N{e) defined 
as follows: 


AJV(e) 


L/2 


lim 

L—> oo 


dx / dy Ap e (x,y ). 


-L/2 


(38) 


FIG. 5. (Color online) 3D plots of electron density distribu¬ 
tion in 2D TI with the defect located at the distance 6 = 3 
from the boundary for two states with different energies close 
to the resonances of the first and second types: (a) e = 0.7 
(first-type resonance) and (b) e = —0.955 (second-type reso¬ 
nance). Other parameters: v = 3.6, a = 2 1 ' 2 , A = 10. 


the other resonance £ e2 was shallow (see Fig. [3]). For this 
potential, we used two energies, one of which e = 0.7 was 
close to the deep resonance and the other e = —0.955 
was close to the shallow resonance. 

Figure [5] clearly shows that the resonances retain the 
main properties of the corresponding bound states in the 
bulk. In the resonances originating from the bound states 
of the first type, the density has a maximum in the cen¬ 


A N(e) gives the spectral density of excess electrons in 
the states flowing around the defect. 

Direct calculations show that AJV(e) has a maximum 
at the resonant energies. The results of specific calcula¬ 
tions of A N(s) are presented in Fig. [b] for different posi¬ 
tions of the defect at a given amplitude of the potential. 
When the defect is located far from the boundary, the 
shape of the peaks on the curve A N(e) is well described 
by the Breit-Wigner formula. When the distance b is less 
than the characteristic lengths of the evanescent states, 
the shape of the peak substantially changes. Neverthe¬ 
less, the integral of AJV(e) over e within the gap is close 
to unity if the peak lies far from the gap edges. 

This fact allows one to estimate the charge accumu¬ 
lated near the defect. If the Fermi energy is close to the 
resonant energy, the accumulated charge is of the order 
of one electron charge. If one assumes that this charge is 
localized in the area of the radius of about 30 nm (this 
is a typical estimation of the decay length of the in-gap 
states in HdTe quantum wells), the potential created by 
this charge can be of the order of 10 mV, which is com¬ 
parable with the gap energy. 

The electron flow in the edge states is strongly dis¬ 
turbed in the vicinity of the defect. To clarify the struc¬ 
ture of the current field in the presence of a defect one 
needs to have an explicit expression for the particle cur¬ 
rent density j in the quantum states defined by the Hamil¬ 
tonian of the BHZ model with the defect. Since the 
Hamiltonian is block diagonal with respect to the spin, 
it is enough to find the current for one of the spin blocks. 
The current density is determined by the term of a diver¬ 
gence in the continuity equation. For the spin-up elec¬ 
trons the current j| is expressed via the spinor compo- 
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FIG. 6. (Color online) Spectral density of excess electrons 
AJV(e) as a function of the energy for a variety of b =3, 3.5, 
4 in two energy ranges: (a) near the resonance of the sec¬ 
ond type, and (b) near the resonance of the first type. The 
parameters used in the calculations: a = 2 1 ' 2 , v = 5, and 
A=10. 


nents %j)\{x,y) and ^{x^y) as follows: 


2 i 


ZD n\ i* 9 ^ 1 (n n\ i* 9 ^ 2 
( B + D X: t — - ( B - D »- 2 — 

o a 

+ —Re [ipiipz] > (39) 


Jfv ~ 


= - 2< Im 


, o , n \ , * (n i~i\ i * 9l ^ 2 

2 a 

+ —Im [ipiipz] . 


(40) 


When the defect is located in the bulk, a circular elec¬ 
tron current is present in each bound state with a given 
spin^S. Its direction is locked to the spin as in the edge 
states. However, there is an essential difference from 
the edge states. The edge state can be occupied by two 
counter-moving electrons with opposite spins so that the 
total current of the filled state is zero. In contrast, in 
the case of a point defect, the problem of the electronic 
structure of the bound state with two electrons, and even 
its very existence, requires a separate study taking into 
account the electron-electron interaction. This issue is 
beyond the scope of this paper. Within the present ap¬ 
proach we study the one-electron states. If one makes 
the natural assumption that the two-electron state has a 
higher energy than the one-electron state, then the one- 
electron bound state is realized. In this state, there is the 
electron current, whose direction depends on the spin of 
the trapped electron. 

If the defect is coupled with the boundary, the field 
of the current density in the edge state flowing around 
the defect includes both the current circulating around 
the defect and the edge current. The configuration of the 
current field is calculated using Eqs. (39) and (40). The 
results are presented in Fig. [7] for two states with the en¬ 
ergies lying near the resonances of the first and second 
types. In both cases there is a current circulating around 
the defect and an edge current flowing around it. As 
the defect approaches the boundary, the circulating cur¬ 
rent decreases and the current flowing around the defect 
increases. 
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FIG. 7. (Color online) Vector plot of the electron current flow 
of spin-up electrons in the presence of the defect located at the 
distance 6 = 3 from the boundary for two states with energy 
(a) e=0.7 near the resonance of the first type and (b) e=-0.955 
near the resonance of the second type. Other parameters are 
v = 3.6, a = 2 1/2 , A=10. 


VI. SUMMARY AND CONCLUDING 
REMARKS 

We have studied the in-gap electronic states induced 
by a nonmagnetic defect with short-range potential in 2D 
TIs and trace their evolution as the distance between the 
defect and the boundary changes. 

If the defect is located far from the boundary in the 
bulk, there are two bound states localized at the defect. 
They exist for both positive and negative potentials of 
the defect. The states are classified as electron-like and 
hole-like states depending on the pseudospin (orbital) 
component of the wave function which vanishes at the 
center: in the electron-like states ip 1 7 ^ 0 and ip 2 = 0 
while in hole-like states i/j 1 = 0 and i/j 2 7 ^ 0 - I n their 
turn these states are classified into two types depending 
on the spatial distribution of the particle density and the 
densities of the pseudospin components. In the states of 
the first type the particle density has a maximum in the 
center while the second-type states are characterized by 
the presence of a minimum of the density in the center. 
In the states of both types there is a particle current cir¬ 
culating around the defect. Its direction is locked to the 
electron spin. 

The presence of two bound states of different types 
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at a given defect potential is a feature of 2D TIs. In 
topologically trivial insulators, there is only one state. 
The existence of two states raises an interesting question 
about the two-particle bound state. 

Another interesting property of the bound states in the 
bulk of 2D TIs is the singular dependence of their energy 
on the defect potential amplitude |u|. The energies of the 
electronlike and holelike states tend to the corresponding 
limiting values e e and as |i>| —> oo. This fact could lead 
to a non-trivial consequence in the case where the crys¬ 
tal contains many different defects with potentials that 
are scattered over a wide range. Since the energy of the 
strong defects slowly changes with their potential, the 
bound state energies are concentrated in narrow spec¬ 
tral bands near e e and £h- The states may overlap and 
form a hopping system, which could manifest itself in the 
transport. 

When the defect is coupled with the boundary, the 
edge states of the host crystal and the bound states trans¬ 
form into a unique set of the edge states flowing around 
the defect. These states have two resonances correspond¬ 
ing to the two types of the bound states in the bulk. The 
resonances retain distinctive properties of the bulk bound 
states. Particularly, they are classified as the resonances 
of the first and second types. The states flowing around 
the defect with the energy close to the first-type reso¬ 
nance have a maximum of the electron density at the 
defect location point, while in the states with the en¬ 
ergy near the second-type resonance the density reaches 
a minimum at the defect. 


When electrons occupy the edge states flowing around 
the defect, the excess charge is accumulated in the vicin¬ 
ity of the defect. Its magnitude can be as large as one 
electron charge. The potential created by such a charge 
can be comparable with the energy gap in 2D TIs. This 
estimation shows that the electron-electron interaction 
can essentially modify the defect-induced states. 

In the edge states flowing around the defect there is an 
electron current for each spin orientation. The current 
density field has two components one of which circulates 
near the defect and another flows around it. The total 
current flowing around the defect depends on the filling 
of the states. In the case of a narrow resonance, even 
a small difference in the population of the states with 
opposite spins can lead to a noticeable current in the 
loop around the defect and the formation of a magnetic 
moment. Estimations show that the magnetic moment 
can be as large as one Bohr magneton. 

An interesting feature of the bulk bound states is the 
presence of a point where the bound states of different 
types are degenerate. This happens at a certain potential 
of the defect. The degeneracy of the resonances is lifted 
due to the interaction of the defect with the boundary. 
The interference of the resonances results in the forma¬ 
tion of a bound state embedded into the continuum of 
the edge states flowing around the defect. 
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